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^ ; ABSTRACT 

Oh' In the framework of the Magnetism in Massive Stars (MiMeS) project, a HARPSpol Large Program at the 3.6m-ESO telescope has 

recently started to collect high-resolution spectropolarimetric data of a large number of Southern massive OB stars in the field of the 
Galaxy and in many young clusters and associations. In this Letter, we report on the first discoveries of magnetic fields in two massive 
stars with HARPSpol - HD 130807 and HD 122451, and confirm the presence of a magnetic field at the surface of HD 105382 that 
^ ' was previously observed with a low spectral resolution device. The longitudinal magnetic field measurements are strongly varying for 

^ ' HD 130807 from —100 G to -700 G. Those ofHD 122451 andHD 105382 are less variable with values ranging from —40 to -80 G, 

and from — 300 to -600 G, respectively. The discovery and confirmation of three new magnetic massive stars, including at least two 
He-weak stars, is an important contribution to one of the MiMeS objectives: the understanding of origin of magnetic fields in massive 
^ ' stars and their impacts on stellar structure and evolution. 

Key words. Stars: massive - Stars: magnetic field - Stars: chemically peculiar - Stars: individual: HD 122451, HD 105382, 
HD 130807 
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1 . Introduction netic fields with t he stellar structure, environment and evo lution 

at high mass (e.g.' Grunhut et al.ll2009l:IOksala et alJl201 ih . 

'""^ ■ MiMeSH (Magnetism in Massive Stars) is a large collaboration onm »u i ■ . uAiDnc i • • j * 

— — I . , . , , . , . In 2010, the polarimeter HARPSpol was commissioned at 

. that aims to address many issues concerning the magnetism of 3.6m-ESO telescope (La Silla, Chile). For the first time, we 
> massive stars. One goal in particular is to determine the global ^^^^^^ ^^^^^^^^ hemisphere with data quality similar 

I > magnetic properties of massive stars with the help of Large . r>i/-v c j xr i -ru * t n 

S/^ /T . . , . . . • . rr- • to ESPaDOnS and Narval. Therefore, a Large Program was es- 

rN Programs (LP) that have been allocated on the high-efhciency . ur u j » i » »u ccn i->/^ c/xr i c u i j 

L_i , . , , . , . _ _ „ , „ ^ tablished to complete the ESPaDOnS/Narval held sample, and 

^ high-resolution spectropolarimeters ESPaDOnS (Canada France ^^^^ ^j^^ ^^.^^ ^^^^^ observing massive stai's in various 

■ ■ Hawaii Telescope, Hawaii) and Narval (Telescope Bernard Lyot, , . , ... t j-a » » • » 

„ s r„, . , , . v^T, open clusters and associations of different ages, to investigate 

France). These programs aim to observe about 200 massive OB .,: u w j »u • » e u 

„ , , , , ^ ^ o^s- , the magnetic held evolution, and the impact of magnetic nelds 

field stars (the Survey Component or SC), m order to search for stellar evolution 
magnetic fields, confirm th ose previously suspected, and deriv e 

statistical properties (Wadi^LaD l2009t iGrunhut et al.l I20TTI) . The HARPSpol sample is divided in two components 

They also aim to observe intensely about 30 aheady known mag- (SC and TC) to follow the same strategy as the Nai'val and 

netic massive stars (the Targeted Component or TC) in order to ESPaDOnS LPs. The HARPSpol SC sample contains about 180 

map in detail their surface magnetic fields. This smaller sam- stars including -110 stars in 7 clusters, and -70 stars in the field 

pie of stars is dedicated to the study of the interplay of mag- of *e Galaxy. The former have been selected from the Catalogue 

of Open Cluster Data (Kha rchenko et al J 12005). while the lat- 
ter have been chosen from the International Ultraviolet Explorer 

* Based on observations collected at the European Southern (lUE) data archive but also from other catalogues or publica- 

Observatory, Chile (Program ID 187.D-0917) tions containing highly probable magnetic stars, in accordanc e 

' http://www.physics.queensu.ca/~wade/mimes with the ESPaDOnS/Narval target selection dWade et alJl2009h . 
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Table 1. Log of observations. Columns 1 and 2 give the date, 
Universal Time (UT) and Heliocentric Julian Date of the obser- 
vations. Columns 3 and 4 give the total exposure time and the 
number of polarimetric sequences. Columns 5 and 6 give the 
peak S/N per CCD pixel (at -501 nm for HD 130807 and at 
-518 nm for HD 122451 and HD 105382) in the specti-a, and 
the S/N per 1.4 km.s ' (for HD 130807 and HD 105382) and 
4.2 km.s ' (for HD 122451) pixels in the LSD Stokes V pro- 
files. Columns 7, 8 and 9 give the longitudinal magnetic field, 
the magnetic detection probability, and the detection type (see 
text). 



Date (d/m) 


HJD ?exp # 


S/N S/N 


Be 


-Pdct 


UT 


(2 455 000+) (s) 


(LSD) 


(G) 









HD 130807 








23/05 05:23 


704.724 


4000 1 500 4600 


2924 


:26 


1.00000 DD 


27/05 06:06 


708.754 


1200 1 520 4800 


-94±26 


1.00000 DD 


28/05 05:58 


709.748 


6000 2 630 5700 


677±21 


1.00000 DD 






HD 122451 








23/05 04:16 


704.677 


1340 5 1680 33200 


-43d 


:20° 


1.00000 DD 


27/05 05:21 


708.723 


1200 20 3000 61800 


-83d 


:14° 


1.00000 DD 


28/05 04:34 


709.690 


1680 14 1260 26000 


-66d 


:29° 


0.98173 ND 






HD 105382 








25/05 01:19 


706.555 


3200 1 990 9650 


-622d 


:26 


1.00000 DD 


26/05 01:59 


707.582 


4000 1 900 8840 


-298d 


:32 


1.00000 DD 


28/05 00:10 


709.506 


4800 1 920 8900 


-406d 


:32 


1.00000 DD 



Notes. Those values are for HD 122451 B only. 



This HARPSpol Large Program was allocated four separate 
runs over two years. During the first run in May 201 1, 57 stars 
were observed including one magnetic calibrator and one TC 
target for which the results will be presented in a forthcoming 
paper. In this letter, we report on the first discoveries of mag- 
netic fields in massive stars with HARPSpol in HD 130807 and 
HD 122451, and we confirm the magnetic field in HD 105382 
previously detected with the low-resolution spectropolarime- 
ter FORS 1 (Kochukhov & Bagnulo 2006; Hubrig et al. 2006). 
Among the 55 SC stars observed during this run, those three 
stars are the only ones in which a magnetic field was detected. In 
Section 2, we present the observations and reduction techniques. 
In Section 3, we detail the HARPSpol results on each star, and 
discuss them in Section 4. 



2. Observations 



We used the HARPSpol polarimeter (iPiskunov et alJ 
combined with the HARPS spectrograph (iMavor et al.l 



2011 



2003) 



installed at the 3.6m ESO telescope at La Silla Observatory 
(Chile), yielding spectra with resolving power /i/A/l of about 
105 000, and covering the 380-690 nm wavelength region. 
All spectra were recorded as sequences of 4 individual sub- 
exposures taken in dififerent configurations of the polarimeter, in 
order to yield a full circular polarisation analysis, as described 
by Donati et al. (1997). The d ata were reduced us ing the pack- 
age "REDUCE" described bv lPiskunov & Valentil (l^2l) . After 
reduction, we obtained the intensity Stokes / and the circular 
polarisation Stokes V spectra of the stars, both normalised to the 
continuum. A null spectrum (A^) was also computed in order to 
diagnose spurious polarisation signatures, and to help to verify 
that the signatures in the Stokes V spectrum are of stellar origin. 
The log of the observations is presented in Table 1 . 



To increase the effective signal to noise ratio (S/N) of 
our data, we applied the Least Squares Deconvolution (LSD; 
Donati et al. 1997) procedure using tailored line masks of ap- 
propriate temperature and gravity for each star. The masks were 
first c omputed using Kurucz ATLAS 9 models of solar abun- 
dance (lKuruczlll993h . with intrinsic line depths larger than 0. 1 . 
We then excluded from these masks hydrogen Balmer lines, and 
lines whose Lande factor is unknown. Finally we have modified 
the line depths to take into account the relative depth of the lines 
of the observed spectra. The resulting masks contain 394, 592, 
and 394 lines for HD 1030807, HD 122451 and HD 105382, 
respectively. The S/N of the LSD Stokes V profiles is about 10 
times larger than the S/N in the original spectra (Table[TJ. 

In order to perform a reliable magnetic field diagnosis, we 
have computed the detection probability inside the LSD V pro- 
files (as described in Donati et al. 1997). We consider that an 
observation displays a "definite detection" (DD) of Stokes V 
Zeeman signature if the probability is larger than 0.99999, a 
"marginal detection" (MD) if it falls between 0.999 and 0.99999, 
and a "null detection" (ND) otherwise (see Table 1). All obser- 
vations of HD 130807 and HD 105382 display DD while two 
DD and one ND have been obtained for HD 12245 1 . The LSD /, 
y, and profiles are plotted in Fig.[T] In almost all of our obser- 
vations Zeeman signatures, as broad as the / profiles, are clearly 
detected in the V profiles, while the profiles are consistent 
with the noise. These results allow us to confidently affirm that 
magnetic fields are present at the surface of these stars. 

We measured for each observation the line-of-sight compo- 
nent of the magnetic field averaged over the visible stellar sur- 
face (the so-called longitudinal magnetic field or B[), by inte- 
grating the / and V profiles over the ranges [-50, 60], [-80, 100], 
and [-70, 105] km.s 'for HD 13080 7, HD 122451, and HD 
105382, respectively (as described by lAlecian et al ] l2009h . The 
values are reported in Table 1 . 



3. Results 

3.1. HD 130807 

HD 130807 (o Lup) is member of the Sco-Cen association 
(iKharchenko et alJ2005h . A companion was detecte d at an angu- 
lar distance varying from 0.07 to 0.14 arcsec ( Perrvman & ES aI 
1997; McAlisteret al. 1990). According to the angular separa- 
tion, the light of both components entered the HARPSpol fibers 
during our observations of this target. 

From a visual inspection we find that most of the spectrum 
of HD 130807 is consistent with a synthetic spectrum of a sin- 
gle star of effective temperature Tetf = 18000 K, surface gravity 
logg = 4.25 (cgs), broadened by vsin/ = 25 km.s"', calculated 
using TLUSTY non-LTE atmosphere models and the SYNSPEC 
code (lHubenvlll988L iHubenv & Lanzlll992h . However, we ob- 
serve that all He i lines are substantially weaker than the syn- 
thetic ones calculated with solar abundance (Fig. |2]i, while the 
Si 11 lines are considerably stronger. The Si, N and Fe lines show 
variability in depth and shape on a timescale of 1 d. These char- 
acteristics suggest that HD 1 30807 is an He- weak sta r with abun- 
dance spots on its surface (iJaschek & JascheO 19741) . 

Magnetic signatures are detected in almost all the lines of the 
spectrum, similar to the LSD one (Fig.[T]/e/f). Many additional 
lines are observed in the spectrum that could be due to Fe ii, 
Fe III, or Ti ii enhancements. All these lines show Zeeman sig- 
natures similar to the others, with the same variations from one 
night to the other. They can therefore be attributed to the same 
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Fig.l. LSD V (top), (middle), and / (bottom) profiles of HD 130807 (left), HD 122451 {middle), and HD 105382 {right). The 
mean error bars are plotted next to each profile. The V and A^ profiles have been shifted and amplified for display purpose. The 
dotted vertical lines indicate the integration ranges for the calculation of B[. left: full black line: 28/05, red-dashed line: 27/05, green 
dot-dashed line: 23/05. middle: full black line: 27/05, red-dashed line: 23/05, green dot-dashed line: 28/05. right: full black line: 
25/05, red-dashed fine: 28/05, green dot-dashed line: 26/05. 



star, rather than a companion, and they are probably the result of 
the chemical peculiarities at the surface of the star 

A significant shift in radial velocity (~ 6 km/s) is detected in 
the strongest spectral lines including Balmer lines, between May 
22 and May 26-27. The maximum reported angular separation 
between both visual components implies a distance > 17 AU, 
and therefore a period > 27 years. This radial velocity shift can- 
not therefore be due to the reported visual companion. A third 
companion very close to the primary could explain these varia- 
tions, but more observations are required to fully understand all 
the peculiarities observed in the spectrum. 

The variations observed in the V profiles over 6 days (Fig. 
[U left) can be understood in terms of the oblique rotator (OR) 
mod el that consis ts of an inclined dipole placed inside a rotating 
star (IStibbslll950l) . The rotational modulation of the shape of the 
V profiles and of the Be values that vary from -94 to 677 G (Table 
1) suggest that the rotation period of the star should be between 
1 and 6 days. More observations, well sampled over the rotation 
period, are required in order to fully characterise the magnetic 
field and better constrain the period of HD 130807. 

3.2. HD 122451 

HD 122451 {p Cen) is a double-lined spectroscopic binary with 
components of similar efifective temperatures (25000 K) and 
gravities {\ogg - 3.5, cgs), and a j3 Cep-type pulsating primary. 
The system is hig hly eccentric (e = 0.835) and orbits with a pe - 
riod of 357 davs ( Ausseloos et al.l200il2006HDavis et al.l2005h . 

We obtained three observations of HD 122451 during 3 dif- 
ferent nights. In order to avoid potential false magnetic detec- 
tions due to pulsations, we have split each observation into many 
sequences of four spectra (Table[T]i, so that the exposure time for 
one sequence is much shorter than the pulsation period. 

The spectra of HD 122451 clearly show two components 
of similar temperatures but different broadening, confirming the 
SB2 nature of the sy stem. We adopt the same definition as 
lAusseloos et al.l (l2006h for the primary and secondary, i.e. as the 
broad-line and narrow-line components, respectively. In order to 
measure the v sin ; of both stars, we performed a least-squares 
fit to few individual spectral lines with the sum of two func- 
tions calculated as the convolution of a Gaussian of instrumen- 
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Fig. 2. Spectrum (full black) of HD 130807 plotted around 
He I 4471 A and Mg ii 4481 A. Synthetic spectra of 18000 K 
(dashed red) and 17000 K (dot-dot-dot-dashed blue) are over- 
plotted. 



tal width and a rotation function as described by 'GraV (19921) 
(see details of the fitting procedure used by Alecian et al. 20081). 
We find a vsin/ of 190 + 20 km.s"' and 75 + 15 km.s"' for 
the primary and secondary, respectively. When compared with 
TLUSTY/S YNSPEC synthetic spectra our observations are con- 
sistent with TefF - 2500 K and log^ = 3.5 (cgs), in agreement 
with the work of lAusseloos et al. (2006). The spectral lines ap- 
pear distorted and show rapid variations very likely due to /? Cep- 
type pulsations. No obvious abundance peculiarity, nor manifes- 
tation of circumstellar matter is observed within the spectra. 
In Fig. [T] {middle), we superimposed the LSD /, V, and 
profiles of our obse rvations. According to the ephemeris of 
lAusseloos et al.l (l2006l) . the 3 observations are roughly at the 
same orbital phase (~0.5), and both components have similar 
radial velocities (~14 and ~4 km.s"' for the primary and sec- 
ondary respectively), which explains why it is diflicult to distin- 
guish both components in the profiles. The shape of the LSD / 
profile shows variations during the run that can be understood 
in terms of radial pulsations in the primary, which would occa- 
sionally broaden the profile. As a result, both components can be 
clearly distinguished in the profile of May 27 (full black line in 
Fig.\l\middle), while it is less obvious in the other observations. 
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Zeeman signatures are detected in many individual spectral 
lines, as in the LSD V profiles. The signatures are as broad as 
the secondary profile, meaning that the magnetic field is detected 
only in the secondary component of the system. However, con- 
sidering the faint Zeeman signatures in the secondary, and the 
broad line shape of the primary, a magnetic field of the same 
strength as the secondary's could exist in the primary, without 
being detected in our observations. In order to estimate the B( 
values of the secondary, we need to extract the V and / profiles 
of the secondary only, from the profiles of the binary. Without 
any evidence of a magnetic field in the primary, we neglect its 
contribution to the V profile, which we consider to be entirely 
from the secondary. On the contrary, the / profile needs to be 
corrected. With this aim, we have first fitted the / profile of the 
binary with the method described above for the individual spec- 
tral lines. Then we have subtracted from the observed / profile 
the fit of the primary. Finally we have measured B( using the 
corrected / and the original V profiles. The values are reported 
in Table 1. 



3.3. HD 105382 

HD 105382 HR 4618) is r nember of the Sco-C en associa- 
tion (iKharchenko et al.ll2005h . iBriquet et all (|2004 classified it 
as He-weak with He patches enhanced where Si is depleted, and 
derived a T^g of 17400 + 800 K, a log^ = 4.18 + 0.20 (cgs), a 
rotation period of 1 .295 + 0.001 d and an inclination angle of the 
rotation axis to the line of sight / = 50 + 10°. 

In our three spectra, we observe strong variations in the spec- 
tral lines, mainly in He i. Si ii and Fe m, tha t are due to abundanc e 
spots on the stellar surface described by iBriquet et al.l (l2004h . 
Clear Zeeman signatures are detected in the metallic and Balmer 
lines, as well as in the LSD V profiles {Fig.\l\righf). The rota- 
tion phases of our observations, calculated with a rotation period 
of 1.295 d, are very different (0.35, 0.14, and 0.63), and yet the 
V profiles are all similarly negative (Table 1). According to the 
OR model, this implies that the magnetic obliquity angle (with 
respect to the rotation axis) cannot be very high (|/S| < 40° if 
i = 50°), otherwise the positive magnetic pole would sometimes 
appear on the visible stellar hemisphere, creating a positive pro- 
file at least once during the run. 

HD 105382 was independently dis covered as magneti c 
by iKoch ukhov & Bagnulo (200"3) and iHubrig et all (l2006h . 
IBriquet et al.i (2007) derived the longitudinal field from FORS 1 
observations and found values ranging from -923 G to 840 G. 
Among their four values, the May 2004 observation (840+58 G) 
is clearly inconsistent with our data as positive values are not ex- 
pected. Pagnulo et al. (2011) very carefully re-reduced the same 
FORS 1 data and found values consistent with those of Briquet 
et al. (2007) except for that observation (for which they derived 
Be = -29 + 69 G). 

We performed a least-square sinusoidal fit to our B( values 
simultaneously with the Briquet et al. (2007) data and could 
find a solution only by removing the May 2004 datapoint. We 
also performe d an in dependent fit using the re-reduced data of 
iBagnulo et al ] d2011h and found a similar result. In both cases, 
the derived period is consistent with that of Briquet et al. (2004). 
The fitted Br values are very similar in both cases, varying 
from -670 G to -20 G, and implying a magnetic obliquity of 
~ 38 ° and a polar field str ength of ~ 2.3 kG, assuming a dipole 
field dBorra & Landstreetl[T980l) . 



4. Discussion 

We report direct detections of magnetic fields in three hot B-type 
stars (18000 K - 25000 K), among a sample of 55 stars in which 
we were searching for magnetic fields with HARPSpol. Two of 
them (HD 122451 and HD 130807) are completely new detec- 
tions. For the other one - HD 105382 - this is the first direct de- 
tection of a Zeeman signature. One of the main MiMeS results is 
the systematic detection of chemical peculiarities at the surface 
of magnetic hot stars (and conversely, e.g. Grunhut et al. 201 lal). 
Among the three stars discussed in this paper two are unambigu- 
ously He- weak. The third one (HD 122451) belongs to a binary 
system with a /? Cep primary, that makes the interpretation of the 
spectrum and the detection of peculiarities inside spectral lines 
very difficult. More observations well sampled over the orbital 
period of the system are required in order to first confirm a mag- 
netic detection in only the secondary, and then disentangle the 
pulsation and chemical peculiarity effects. 

The interplay between radiative and magnetic forces 
([Hunger & Groote 199 9) is usually assumed to be at the origin of 
the over- or under-abundant He spots at the surface of hot mag- 
netic B stars. These spots are very often correlated with the stel- 
lar magnetic f ields (e.g. Veto 1990). In the case of HD 105382, 
IBriquet et all ( 12004) found a large He spot at a latitude of 60°. 
If our estimate of its magnetic obliquity is confirmed, this spot 
would be situated close to the South magnetic pole, demonstrat- 
ing once more the importance of magnetic fields in the formation 
of chemical spots. 

More observations of these magnetic stars are planned within 
the HARPSpol large program in order to perform detailed map- 
ping of their magnetic fields, and better confront the models and 
theories of magnetic massive stars. 
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